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Recent reports of excess high-Q^ events in both the HI |I[] and ZEUS detectors at 
HERA are of considerable interest, since they correspond to possible leptoquark sightings. 
The HI Collaboration |jl| has reported 12 neutral current (NC) events with Q^ > 15000 
GeV^ for an integrated luminosity of 14.2 pb^^ against a Standard Model (SM) deep 
inelastic scattering (DIS) background of 4. 7 ±0.76 events. The fluctuation probability for 
this is 9 X 10-3. The ZEUS Collaboration has reported 5 NC events with Q^ > 20000 
GeV^ for an integrated luminosity of 20.1 pb^^ against a Standard Model background 
of 0.9 ± 0.08 events, for which there is a fluctuation probability of 1.4 x 10^^. The 
excess events appear to indicate a resonance of mass around 200 GeV. A confirmation 
would point emphatically to new physics since leptoquark fields — vector or scalar — 
are not present in the SM. Vector leptoquarks arise naturally in Grand Unified Theories 
(GUTs) 0, but these are constrained by the proton lifetime to be very heavy — far 
beyond the reach of HERA energies. Relatively light vector leptoquarks may still exist in 
certain models, but these are not readily amenable to coupling constant unification |Q. 
On the other hand, scalar leptoquarks — arising from Higgs multiplets — are generically 
free from such constraints. Most such models, however, suffer from a major conceptual 
drawback, namely the lack of a satisfactory explanation for the large mass hierarchies 
in the theory. The presence of a light leptoquark accessible at HERA merely aggravates 
this problem. In contrast, the best known answer to the problem of mass hierarchies 
is provided by supersymmetry p[, which is also known to facilitate coupling constant 
unification p. It seems, therefore, that one will eventually require supersymmetry as 
well as grand unification in order to have a satisfactory theory involving the presence of 
elementary leptoquarks. However, before looking for an explanation of the HERA events 
within the framework of any supersymmetric GUT, it is interesting to ask if these events 
can be accommodated within the minimal supersymmetric extension of the Standard 
Model (MSSM) ^, which is a simpler construction. 

In its usual formulation, the MSSM respects the global conservation laws for the 
baryon number B and lepton number L, which are valid in the SM. More conventionally 
represented as conservation of -R-parity, Rp = (_ i)3B+I'+2S'^ where S is the intrinsic 
spin [0, this feature immediately precludes leptoquark-like interactions. However, this 
symmetry was originally imposed ad hoc and has no compelling theoretical motivation. It 
is thus important to investigate the consequences of possible i?p-violating terms |^, |^ in 
the Lagrangian. This assumes particular significance in view of the fact that it would allow 
the squarks to have leptoquark-like interactions. A broken i?-parity (^p) would also imply 
that the lightest supersymmetric particle (LSP) is no longer stable: search strategies and 
existing bounds on MSSM masses and couplings need to be modified accordingly. 

While it is possible to have both L- violating and 5- violating ^pterms, non-observation 



of proton decay imposes very stringent conditions on their simultaneous presence |]T0 
Assuming that the baryon-number- violating terms are identically zero 0] helps evade this 
constraint in a natural way, apart from rendering simpler the cosmological requirement 
of the survival of GUT baryogenesis through to the present day. This survival can then 
be assured if at least one of the lepton numbers L, is conserved over cosmological time 
scales O . 



We thus concentrate on the L-violating part of the superpotential, which can be 
written (with the MSSM superfields in usual notation) as 



"^ ~ T^^ijkLiLjEk + X^^i^LiQjDk, 



(1) 



where Xijk = —Xjik-, with i, j, k being family indices. Since we are interested in processes 
involving quarks, we shall further restrict ourselves to the A'-type couplings only. The 
absence of Xijk couplings can be assured by the imposition of suitable discrete symmetries 



ri| . Neglecting quark mixing effects, the interaction Lagrangian can then be expressed 



as 



-A' 



ijk 



T^ihdkRdjL + djidkRi^iL + dljiii^iiYdji 



-diidkRUjL - UjidkRCiL - dl^^enYuji + h.c 



(2) 



Clearly, while the sleptons and sneutrinos behave rather like charged and neutral Higgses 
in a non-minimal model, the squark coupling mimics that of leptoquarks. 



Given the above couplings, the processes of interest at HERA are ||13|, |1^, |T^ 
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(3b) 

and, of the 27 possible A's, only nine (A'^^^-^) are relevant to our discussion. If one admits 
the possibility that more than one of these couplings may be non-zero, tree-level flavour- 
changing NCs are introduced and the bounds on most of the products are very severe IJTB 



More conservative bounds are obtained under the simplifying assumption that only one 
of these couplings may be non-zero. The relevant 2a bounds (with their sources) are : 
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[la) Atomic parity violation 117] 



Now, of the nine couplings X[jk, only for three (A'^^^^^) can a valence quark in the proton 
participate in the production mechanism (eq.(^)). For all other couplings, a sea-quark 
must appear in the initial state. As the data ||l], 0] seem to indicate a moderately large 
value for the parton momentum fraction, the corresponding sea-quark densities are rather 
small. An explanation of the data would then need a large value of the corresponding A', 
which is already ruled out by the constraints of eq.(|l). 



Of the three that are left, X'm has aheady been constrained to be too small to be 
relevant. Consequently, we are left with two possible couplings and hence scenarios, 
namely : 

A'i2i • e"*" + (i — > cl — > e~^ + d (5a) 

and 

A'lgi : e+ + rf — >iL — > e+ + d, (5b) 

so that the HERA NC events could be signals for the production of a left-handed scalar 
charm or stop. The major difference between these signals and those for a leptoquark (of 
the same mass and quantum numbers) lies in the fact that these squarks can have other 
(i?p-conserving) decay modes into quarks and charginos/neutralinos. This would show up 
either through the branching ratio or the width of the resonance ||22||. 



Until now, we have considered only the ^3 = | component of the squark doublet. How- 
ever, the masses of the two isospin components are related, with the splitting determined 
by the corresponding quark masses and the S'?7(2)i-breaking D-teim 0: 

2 2 2 I 2 2 na 

rrij^ = m^^ — m^ + mf, — my^ cos ip , 

and similarly for the (c/,, sl) pair. This implies the presence of at least one more relatively 
low-lying squark state {sl is somewhat heavier than cl while hi is somewhat lighter than 
II)- a possible consequence is an extra contribution ||2^ to 5p = 1 — m^/m| cos^ 9w- It is 



easy to check that the contribution due to {cl, sl) is much smaller than the experimental 
errors |2^. For the stop solution, though, this constraint is non-trivial, especially for 
lighter stop masses within the (180-220) GeV range. As an example, for m^^ = 180 GeV 
and rrit = 170 GeV, one needsQ tan /5 ^ 4 for consistency with 5p at the la level. This also 
implies that the 6l would be just beyond the energy range of LEP2. 

As far as HERA is concerned, the two processes (^ & |5E| ) are almost identical, the 
only major difference being the fact that the neutralino decay modes of the stop are 
greatly suppressed by phase-space considerations. Consequently, for most of our analysis, 
we shall consider ( |5a| ) to be the representative process and will highlight differences only 
where they are crucial. 

While the production cross section for squarks is obviously quadratic in A', the sub- 
sequent decay can be more comphcated. If i?-parity is violated, the squark may indeed 
be the LSP without contradicting any phenomenological bounds. In this case, it would 
be virtually indistinguishable from a leptoquark coupling to left-handed electrons and 
the right-handed (i-quark. However, we do not make any such assumption and allow the 
squark {q^ = cl or ti) to have i?p-conserving decays : 

Ql — ^ g + X? , 

(6) 

Ql — > q' + xt ■ 



^This applies for the case of vanishing left-right squark mixing. In the presence of mixing, the con- 
straint may be considerably relaxed. 



The various branching fractions depend on both the size of A' and the quark-squark- 
gaugino couphng as well as on the gaugino masses. Under the assumption of gaugino 
mass unification in the electroweak sector, the latter are determined in terms of three 
parameters: the gaugino mass term M2 in the SU{2)i sector, the higgsino mixing param- 
eter /i and tan (3, the ratio of the Higgs vacuum expectation values. Any explanation of 
the HERA events within the framework of -Rp- violating supersymmetry must, therefore, 
consider the constraints that can already be imposed on the parameter space from LEP 
data. Now, if i?-parity is violated, most of the usual bounds do not hold any longer unless 
further assumptions are made. In fact, with our choice of a non-zero A', the only relevant 
constraint (in the M2-/i plane) is the one derived from the total width of the Z ||25|| : 
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Figure 1: The range in A' that can account for the HERA events as a function of the 
squark mass. The minimum value is common to the two pairs {cl, X121) and (tL, ^'12,1) ■ 
The spread accounts for all values of (M2, fi,ta.n(3) consistent with the LEP constraint 



We make a simple parton-level Monte Carlo simulation of the process e^p -^ ci/ii + 
X — > e^d + X aX HERA energies using acceptance cuts and selection criteria closely 
modelled on those adopted for the NC signal by the HI Collaboration [jl|. Like them, we 
use the MRS(H) [^ distributions for the parton densities inside the proton, although we 
have checked that our results are not very sensitive to the particular structure functions 
being used. For a given value of the squark mass, we determine all values of the ^pcoupling 
X'lji {j = 2, 3), which can yield the observed excess of 12 — (4.7 ± 0.76) ~ 7 events in 14.2 
pb~^ of data |l|. A value of A^^^ is deemed acceptable if we can find a set of values of 
M2,/i, tan/3 in the ranges (0 < M2 < 1 TeV), (-1 TeV < /x < 1 TeV), (1 < tan/5 < 50), 
which is consistent with the LEPl bounds and for which the required number of events is 
predicted. Our results are shown in Fig. |I| where solid (dotted) lines show, for each value 



of the squark mass, the upper and lower permissible values of A'^^^ for j = 2 (3). The lower 
bounds are the same for the two cases; this corresponds to the fact that a small coupling 
requires a large branching ratio into e"*" + d, which is only possible if both charginos and 
neutralinos are very heavy, irrespective of the mass of the quark in the decay products. 
Another consequence of this fact is that the dependence of the allowed range in coupling 
is minimally dependent on tan/3. Thus, even if low values of tan/3 are disallowed by the 
p-parameter constraint for stops, the curves remain unchanged. It is interesting that the 
coupling cannot be smaller that 0.03 or larger than 0.26 (0.22) for cl (t^) production to 
be the explanation of the HERA events. This result is consistent with the earlier bounds 
given by the HI Collaboration |15] based on 2.8 pb^^ of data. 
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Figure 2: The Rp-violating branching fraction for the squark as a function of the coupling 
for three different values of the squark mass (marked, in GeV). 



In Fig. ^ we exhibit the branching fraction of UiL ^ e+ + d (z = 2, 3) required to 
explain the HI events for Q^ > 15000 GeV^ as a function of the relevant A', for three 
representative values of the squark mass. These curves tell us that for reasonably large 
values of the coupling the branching ratio must be rather small, i. e. 0.25 or less (unless the 
squark mass is at the lower end of the range). This immediately points to the conclusion 
that there should be a large cross-section for the squark in question to decay to charginos 
or neutralinos, which, in turn, should be expected to have decays leading to recognizable 
signals. Thus, under the assumption that only one i?-parity-violating coupling is non- 
zero, the lightest neutralino will decay into two quarks and an electron/neutrino. For 
example. 
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121 
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as well as the conjugate modes. The situation is analogous for A'^g^, except for the fact 
that the decay proceeds only through the dxid and bxi^ couplings. Chargino decay is 
more complicated. Apart from the i?p-violating decay (isospin analogue of the process in 
eq.(||)), it could decay into two SM fermions and Xi with the last-mentioned then decaying 
as above. The ratio of the i?p-conserving and i?p-violating decay modes is determined by 



both the size of A' and the gaugino content of the chargino, and the chargino-neutrahno 
mass sphtting. These issues have been studied in detail in ref. |]15[, for an integrated 
luminosity of 2.8 pb^^ with the conclusion that such cascade decays of a squark in the 
(180-220) GeV range cannot be recognized above the background for A ^ V^vra (~ 0.31 
at HERA energies) . Since this is precisely the region of the parameter space that we are 
interested in, we may safely conclude that published constraints from HERA [^ do not 
rule out any part of the parameter space shown in Fig. |l|. On the other hand, with the 
higher luminosity accumulated since then, it may well be possible to observe such signals. 
If, indeed, such signals are seen, the signals will be able to clearly distinguish i?p-violating 
signals from those of leptoquarks. Non-observation would not rule out the ^p option, 
but would restrict us to the lowest range of A' for which the decay mode is practically 
indistinguishable from that of a leptoquark. 
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Figure 3: The number of high Q^ events due to squark production as a function of the 
minimum Q^ used to bin the e"*", for three representative squark masses (marked, in GeV). 
The HI data /0/ are also shown. 

The hypothesis that the excess in high-Q^ events is due to resonant production of 
a particle with leptoquark-like interactions is subjected to a simple test by the data 
presented by the HI Collaboration for different cuts on the minimum Q^. They 
observe an excess over DIS predictions of approximately 7, 4 and 2 events, for Q^ > 
15000, 20000, 30000 GeV^, respectively. This is also consistent with the excess observed by 
the ZEUS Collaboration |^ , given their slightly higher luminosity. In Fig. |^ we illustrate 
our predictions for squark production in i?p-violating supersymmetry, as a function of 
the cut on minimum Q^, normalized to 7 events, which agrees with the HI excess for 
Q^ > 15000 GeV^. As before, solid, dashed and dotted lines correspond to squark masses 
of 180, 200, 220 GeV. The flatness of the initial part of the curve and its linear decrease 
thereafter are indicative of the relative smallness of the t-channel contribution. This is 
as expected, since it arises from a heavy sea-quark (c or t) and is a non-resonant process. 
It is clear from the figure that the observations seem to favour somewhat lower values of 



the squark mass within the range of interest. However, for such low statistics, any such 
statement must be made with a great degree of caution. 

We have also examined the implications of an i?-parity-violating solution, for the 
charged current (CC) excess seen by the HI Collaboration ^. This is rather hard to 
explain from direct squark decays. As the neutralino always decays within the detector 
for the allowed range of A' couplings, the only possibility is to have the following process 

e+ +P^cl + X ^ c + x^,+X -^ c + ssu{u)+X , (9) 

where the neutralino decay takes place through a virtual sl squark, which has to be close 
in mass (see above) to the cl- Both the lightest and next-to-lightest neutralinos may 
participate in this process. The final state would be purely hadronic (mostly with more 
than one jet) with missing momentum, arising from the neutrino. A parton-level Monte 
Carlo analysis shows that for the greater part of the MSSM parameter space the final 
state neutrino is rather soft. In fact, two thirds or more of the events are lost because 
of the selection cut of 50 GeV on the missing momentum which has been imposed by 
the HI Collaboration. For values of A' consistent with the NC data, the branching ratio 
of the Cl to neutralinos can be as high as 70%. For the same set(s) of parameters, 
the branching ratio for the cssi'{u) mode is around 40%. Convoluting these with the 
production cross-section, we predict up to 3 excess events in the CC sample for Q^ > 15000 
GeV^, depending on the choice of the point in the parameter space. The excess drops 
rather sharply for Q"^ > 20000 GeV^, however, which does not seem to agree well with the 
trend of the HI data. Pushing the comparison further is not very meaningful in view of 
the low statistics. It is also noteworthy that the ZEUS Collaboration [§] has not reported 
any CC events, although their luminosity is somewhat higher. A full investigation must, 
therefore, await more data. Meanwhile we can say with some certainty that a charged 
current excess at the level observed so far can be accommodated quite easily within our 
^p-scenario. 

It may be said, therefore, that the observed excess in high-Q^ events at HERA can be 
simply explained by squark production in i?p-violating supersymmetry. However, unless 
a further search for chargino/neutralino decay modes of these squarks is made and actual 
signals found, it will be difficult to decide whether the present signals are those arising from 
leptoquarks or squarks. We turn, therefore, to other possible consequences of the proposed 
scenario, with special emphasis on those that might distinguish it from the leptoquark 
alternative. We make the conservative assumption that no other supersymmetric particle 
would be discovered directly. 

We start with LEPl results on Rb = T{Z -^ hb)/T{Z —>■ hadrons). This is of relevance 
only to the {ti, ^'m) solution. It might be argued that a relatively light stop could 
contribute significantly to this ratio, especially if the chargino were light too. One must 
remember, however, that we are dealing with ti and the effect is smaller compared to that 
for tR. Explicit computations show that the extra contribution is of the order of 10~^ 
to 10^^ and hence is consistent with the observations [27]. If we assume a degenerate 
(ti,t/j) pair (this does not affect the HERA signal), one can get a contribution as large 



as 10 ^, which is also consistent with observations. The contribution due to the A'^3;^ 



couphng itself is small [^. The situation for Ri is similar [^. For X[^i ~ 0.25 or less 



the excess contribution is always well within the experimental errors at la. Essentially, 
contributions to both Rh and Re are severely suppressed by a squark mass in the 200 
GeV ballpark. Thus precision measurements at LEPl will not help to confirm or deny 
our hypothesis. 

Since these squarks are too heavy to be produced at LEP2, one could only consider 
virtual effects. A relevant process is e~^e~ — > cc. A non-zero A'^gi leads to an additional 
(t-channel) diagram. A study of the angular distribution of the charm jets may rule out 



the high A'j^gi ^^^ ^^ the parameter space [Q. It is only fair to point out, however, that 



this is just the range where further searches at HERA itself may be effective. 

At the Tevatron, there is enough energy to produce a pair of these squarks. Both Rp- 
conserving and /2p-violating processes will contribute. The cascade decay modes (through 
charginos and neutralinos) have been studied in the literature [^ in fair detail. The 
upshot of these is that squark masses and couplings of the magnitudes conceived of here 
would not lead to signals distinguishable from the background. On the other hand, direct 
decay into e~^d (in common with the leptoquark signal) will give rise to two hard jets and 
a dielectron pair. This signal might be of significance I^Of, especially for low values of A'. 
However, in view of the large SM backgrounds, a detailed and separate study seems to be 



called for 31 . 



Virtual effects at the Tevatron are not likely to be a decisive discriminator. A non- 
zero X[^i will enhance the ti production cross-section, but the relatively large errors in the 
cross-section measurement [Q do not allow us to probe X'l^^ ^0.3 unless the selectron is 



light enough to be seen at LEP2 |33]. The analogous process for A'^^ji would result in an 



enhancement of the large-pT dijet rates. An improvement in the angular measurement |]3^ 
of the pure dijet spectrum could, in principle, rule out the large A'^^ji end of the parameter 
space. This, however, is unlikely to substantially improve upon the bounds derivable from 
cc production at LEP2. Moreover, as pointed out already, this range of values of A' can 
probably be probed at HERA itself. Dilepton production, on the other hand, can probe 
both the couplings |^, but would not be able to distinguish between them or from the 



leptoquark alternative. Thus, although they are certainly worth a detailed investigation, 
one cannot expect immediate confirmation or otherwise from virtual effects. 

Another interesting possibility, which arises in models with i?p-violation, is that of 
single squark production at the Tevatron, the process relevant to the present investigation 
being g + d —^ Ujl + e~ , where j = 2, 3. With the squark decaying through the ^p 
interaction, we are led to the signal p + p — >• e"'" + e~ +jet + X. However, as Fig. ^shows, 
choosing a large value of A' to get a sizeable cross-section here is more than offset by the 
HERA requirement of a small branching ratio of the squark to e~^d. As a result, the signal 
is strongly suppressed and is not likely to be distinguishable above the SM background. 
It may still be possible to detect signals of the chargino or neutralino decay modes of the 
squark (as at HERA), especially with added luminosity, such as will be available with the 
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commissioning of the Main Injector. Eventually, perhaps, CDF/DO will be in a position 
to confirm the signal or, at least, rule out parts of the parameter space. 

In summary, then, we have investigated the possibility that the excess high-Q^ events 
observed at HERA are due to the production of squarks of i?-parity- violating supersym- 
metry. We have shown that in view of low-energy constraints, this is acceptable only if 
one of the couplings A'^^gi; -^isi i^ non-zero, so that the produced particle must be a left- 
handed scalar charm or stop. The coupling must lie between 0.03 and 0.26 for the first 
option and 0.03 and 0.22 for the second, depending on the mass of the squark, which, 
according to the data, must lie around 200 GeV with an uncertainty of about 10%. It 
turns out that, for a large range of these couplings, large branching ratios of the produced 
squark into charginos and neutralinos are predicted. Decays of the latter should lead to 
distinctive signals at HERA. On the other hand, small values of the -^p-violating coupling 
would make it extremely difficult to distinguish squark signals from those arising from 
leptoquarks of similar masses and couplings. We also show that confirmation (or other- 
wise) of our hypothesis at the two running high-energy facilities, LEP and the Tevatron, 
is unlikely to be immediately forthcoming, although we might expect some results from 
the Tevatron eventually with higher luminosities. The supersymmetric solution would 
thus remain a valid and exciting hypothesis to explain the HERA events, and we look 
forward to new results in this area in the near future. 

The authors are grateful to J. Kalinowski for drawing our attention to the HERA 
events, to M. Drees for useful information and to C.E.M. Wagner for discussions. The 
parton distributions were calculated using the package PDFLIB 
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